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We study the most suitable procedure to measure the effective temperature in off-equilibrium 
systems. We analyze the stationary current established between an off-equilibrium system and a 
thermometer and the necessary conditions for that current to vanish. We find that the thermometer 
must have a short characteristic time-scale compared to the typical decorrelation time of the glassy 
system to correctly measure the effective temperature. This general conclusion is confirmed analyz- 
ing an ensemble of harmonic oscillators with Monte Carlo dynamics as an illustrative example of a 
solvable model of a glass. We also find that the current defined allows to extend Fourier's law to 
the off-equilibrium regime by consistently defining effective transport coefficients. Our results for 
the oscillator model explain why thermal conductivities between thermalized and frozen degrees of 
freedom in structural glasses are extremely small. 



I. INTRODUCTION 

In the last decades there has been an increasing interest in trying to extend the thermodynamics ideas to systems 
which are far out of equilibrium. The general problem is to quantify non-equilibrium fluctuations within a thermo- 
dynamic approach. Recent studies suggest that the slow dynamics of glassy systems allows us to do such kind of 
generalizations!] . In this direction, it has been proposed a definition of a non-equilibriuna temperature (also called ef- 
fective temperature) based on the violation of the fluctuation-dissipation theorem (FDT)a This temperature verifies 
the relation: 

TeS {t',t)R{t\t) = ^^- t>f, (1) 

where C{t' ', t) is the autocorrelation function and R(t', t) is the corresponding response function. It has been suggested 
that for a large class of glassy systems in the long-time regime (t 1 — > oo) the effective temperature is only a-function 
of the smallest time, i.e T e e(t',t) — ► T c g(t') and also that is a good quantity in the thermodynamic senseErEI. This 
knowledge is mainly theoretical and comes from the analysis of mean-field glassy systems. There is also good numerical 
support to these assumptions for structural glass models^, although from the experimental point-qf view, the question 
is still far from settled (see, for example, recent experimental results on glycerolo and laponitaid). Note that in the 
general framework of thermodynamics, the concept of temperature is associated with local equilibrium, i.e the fact 
that local regions in space admit a thermodynamic description in terms of space-dependent intensive variables such 
as the temperature. 

In this paper we try to understand the most adequate procedure to measure the effective temperature. This question 
has been recently addressed by Cugliandolo, Kurchan and PelitiB and Exartier and Pelitllil who have discussed possible 
procedures to measure the effective temperature. 

In a previous paper we have used the concept of the effective temperature in order to investigate the validity of the 
zero-thermodynamic law in off-equilibrium systemsEj. We have seen that in absence of coupling via the dynamics two 
coupled systems far from equilibrium never thermalize. This poses the question: why different degrees of freedom in 
structural glasses can stay at different effective temperatures forever?. In this paper we shall show that the answer 
comes from the smallness of the conductivity between the two systems in the aging state. 

The paper is organized as follows. In section II we define the heat transfer between two generic coupled systems. 
In section III we analyze the heat transfer between an off-equilibrium system and a thermometer in the framework of 
Linear Response Theory (LRT). In section IV we show an explicit example. In section V we extend Fourier's law to 
the off-equilibrium regime. In section VI there are the conclusions. 
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II. DEFINITION OF THE HEAT TRANSFER 



In systems which are in equilibrium, the concept of a thermometer is related to the zeroth law of thermodynamics. 
If we want to measure the temperature of a system in equilibrium with a thermal bath at temperature T, we take 
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another macroscopic system (called thermometer) in contact with its own thermal bath. Then we put the thermometer 
in contact with the system in such a way that the system is weakly perturbed. If the temperature of the thermometer's 
bath is different from T , the whole system (thermometer + system) is off-equilibrium, as a consequence there is a 
heat transfer from the hot system to the cold one. In fact, in order to measure T we change the temperature of the 
thermometer's thermal bath until we cannot measure any heat transfer between the system and the thermometer; 
then, the whole system is in equilibrium and the temperature of thermometer's bath is just T . 

For the non-equilibrium case we shall follow the same steps. We couple the off-equilibrium system with a thermometer 
which is in equilibrium with a thermal bath at a certain temperature. If perturbed, this thermometer relaxes to 
equilibrium with a characteristic time-scale (in this case the perturbation is the coupling with the off-equilibrium 
system). In order to have some kind of "off-equilibrium zeroth law" we have to define the heat transfer between the 
two systems. Let us suppose that we have two generic systems, X and Y with N internal degrees of freedom Xi and 
yi respectively and we couple them. Assuming that the coupling between the systems is small, we may write for the 
total energy: 

N 

H^H-eJ^ZiVi ■ (2) 

i=l 

We consider the two-time correlation and overlap functions defined in the following way: 

N , N 



Cx(M0 = ^I>(*)zi(O; ^(M0 = ^ !>(«') (3) 
i=i i=i 

i N i N 

Qi(*'0 = ^X>i(*)zi(*'); Q2(t,t') = -J^^Mt') ■ ( 4 ) 
i=i i=i 

FollowingB'til we define the heat transfer between them as the net balance of the power supplied by one system to the 
other. In terms of the overlap functions the heat transfer J{t,t') reads: 

8 N 
J(i,f') = e-(Q 2 (t,t') - Qi(f,t)) = e]T (^(t)^') - • (5) 

i—l 

The heat transfer vanishes in equilibrium because of the time-reversal symmetry property of correlation functions, i.e 
Qi(t, t') — Q2(t, t') — Q{\t — t'\). If there is a system which is out of equilibrium this quantity is, in general, non-zero. 
From now on we shall assume that the system X is the thermometer, while Y will be the off-equilibrium system 
which is relaxing towards its equilibrium state. Depending on the ratio between the characteristic time-scale of the 
thermometer and the relaxation time-scale of the off-equilibrium system we can identify the different contributions of 
each term in the right hand side (r.h.s.) of (||): 

• If the relaxation time of the thermometer is fast enough, only the second term in the r.h.s. of M) contributes. 
The contribution of the first term is negligible due to the fact that the net focce made by the fast degrees of 
freedom of the thermometer upon the system is zero on average. Note that inQ the definition of heat transfer 
contains only the second term in the r.h.s of (|^); although not explicitly stated there, in this definition it is 
implicitly assumed the fact that the characteristic time-scale of the thermometer is much smaller than the 
relaxing time-scale of the off-equilibrium system. 

• If the thermometer has a relaxation time larger than the characteristic time of the system, the main contribution 
to J(t,t') comes from the first term in the r.h.s of (^). We will see that in this case we don't have a good 
thermometer for measuring effective temperatures. Obviously, when the typical time of the thermometer is of 
the order of the relaxation time of the system, the two terms in the r.h.s of (j^) arc of the same order, and we 
must to take into account both of them. 

It is believed that for measuring the effective temperature the characteristic time-scale of the thermometer should be 
of the same order as the relaxation time-scale of the system. This is not clear at all because the power supplied by 
the thermometer to the system must be negligible in order to measure correctly the effective temperature, so the term 
x(t)y(t') should be smaller enough in comparison to the term x(t')y(t). This condition can be fulfilled only in the case 
of faster thermometers. Moreover, if the time-scale of the thermometer is longer than the relaxation time-scale of the 
system, the two systems are off-equilibrium during the measurement due to the fact that the perturbed thermometer 
has not time enough to relax to its equilibrium state. In fact, the temperature of the thermometer is only well defined 
when it is in equilibrium. 
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III. LINEAR RESPONSE THEORY 



Now suppose that e is small enough so that we can consider the e-coupling term in (H) as a perturbation. Thus the 
effective temperature will be slightly perturbed. Hence, using linear response theory (hereafter referred to as LRT), 
we have: 



-°(t) 



duRx(t,u) yi (uy, yl(t') = yr°(t') 



duRyit' , U)Xi(u) 



(6) 



Here we have considered that the two systems are coupled at a time t = s (the "coupling time"); so, we must integrate 
the responses of the system and the thermometer between two times: the coupling time s and the time we want to 
measure. For the general non-equilibrium case, if we assume that T e g(t',t) only depends on the smallest time, we 
obtain in the limit when the two times involved are much bigger than the coupling time s: 



lim J = e 2 



dC Y {t,u) 
du 



Cx{t,t)(Px ~ Pv{t)) 



I duRx(t,u 

J s 

duCx(t, u 



Py(u) dC Y (t,u) dC Y {t,u) 



P. 



x 



du 



Ot 



d 2 c Y (t,u) ,d 2 c Y (t,u) 

Px — h Py[u) 



du 2 



dudt 



(7) 



where we have used (||) and the generalization of FDT (|l|), being /3 Y (u) = 1/T cff (u) and (3 X = l/Ijf. Is easy to verify 
that the expression for the current in the equilibrium case reduces to (Q) with j3 Y (u) = j3 Y , yielding: 



lim J(t,t') = 2e 2 (/3 x 
t->e 



* dC x (t 
du- 



ll) dC Y {t-u) 



(8) 
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which is some kind of Green-Kubo relationEl. From this expression we realize immediately that in the equilibrium 
case the heat transfer is zero only when the two systems are at the same temperature, and it remains equal to zero 
for all times. At this point we have to explain better the measurement procedure. As the off-equilibrium system is 
relaxing towards its equilibrium state, we cannot find a single temperature thermometer that cancels the heat transfer 
during all the time. This is due to the assumption that the thermometer is a system in equilibrium while the effective 
temperature of the off-equilibrium system depends on time. In order to measure the effective temperature at a time 
t we have to couple the thermometer at a time s. We shall see later that we can take s = t — tit where r is the 
characteristic time-scale of the thermometer and n a large integer number. Physically the result is independent of n 
provided that n » 1 . Obviously, for every time t we are interested in, the measured temperature will be different 
as we expect. It is easy to prove that only when the characteristic time of the thermometer r is much less than the 
relaxing time of the system, i.e t << t , we have: 



otherwise, for the limit r >> t, we have: 



Um (fc-M t ))_„ 

t^oo f3 Y {t) 



t^oo f3 Y (t) 

and the thermometer cannot properly measure the effective temperature. 



(9) 



(10) 



IV. AN EXAMPLE 



As an instructive example we consider.^ .pinrlel of uncoupled harmonic oscillators with parallel Monte Carlo dynamics 
introduced ino and recently revisited ir£3ll3'E3. This is a simple model which contains the essential features of glasses 
such as aging in the correlation and response functions. This model is defined by the Hamiltonian: 

1 N 

H=-Kj2*i ■ (11) 

i=i 

The Monte Carlo dynamics is implemented by small movements: Xi — > x% + ri/yN where 7"$ are random variables 
Gaussian distributed with zero average and variance A 2 . The move is accepted according to the Metropolis algorithm 
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with probability W(AE) which satisfies detailed balance: W(AE) — W(—AE) exp(— (3AE) : where AE is the change 
in the Hamiltonian. 

Here we will concentrate our attention on the dynamics at zero temperature which is known to be glassy. The interest 
of the zero-temperature dynamics relies on the extremely slow decay of the energy which decays like l/log(t) due 
to the entropic barriers generated by the low acceptance rate. For our purpose, we need the correlation function 

defined as C Y (t,t') — j?Ylf ' ViW)yi(t') an d the response function: R Y (t, t') = ^ g^ffi] ) with t > t', where the 

magnetization forJihis model is M(t) — Yli 
The results areO: 

Cy(t, t') = exp (- £ f Y (x)dx^j ; R Y (t, t') = f -^l exp (- J* f Y (x)dx^j 0(t - t') , (12) 

where f Y (t) is a function that contains the whole dynamics and decays like 1/t at zero temperatureO. The effective 
temperature ([!]) can be exactly computed for this model: 

■ < 13 > 

Note that for this particular case, the effective temperature is only a function of the smallest time. This is a 
characteristic feature of this particular model that is believed to hold for structural glass models in the asymptotic 
limit t' — > oo. In this limit ( |l3| ) verifies T e g(t') — ► 2E Y (t') at zero temperature. 

To measure this effective temperature we have to define a general thermometer (i.e a system in equilibrium which 
verifies FDT). In general we can write for the correlation and the response functions of the thermometer: 



C x (t-t') = C x (0)e- (t -'' )/T ; R x (t-t') = ^C x (t-t') . (14) 



Notice that if the characteristic time of the thermometer r is very small (a fast thermometer) the correlation function 
Cx(t — u) only contributes if u « t ; in this case we can write dcy (t,n) _ — dCy g^' u ^ ■ This result is valid in the t — > oo 
limit and comes from the fact that, at very short time-scales, the aging system can be considered as if it were in 
equilibrium (at the corresponding effective temperature) and the partial derivatives of the correlation functions can 
be interchanged in the usual way. This can be checked by explicit derivation of eq. (|l2"|) . Then, the heat transfer 
will vanish if the condition (3 Y (t) — (3x is satisfied. Therefore, if we have a very fast thermometer we can measure 
exactly the effective temperature. By numerical evaluation of the integrals (Q) we compute the temperature T m of 
the thermometer which makes the heat transfer vanish. The results for different thermometers are shown in the left 
panel of Fig.l where one can clearly see that for t » r the measured temperature T m = l/p m is always higher than 
the effective one and converges to it logarithmically. On the contrary, for t « r the thermometer does not measure 
(for a given time t) the right effective temperature. We should say that for t < < t the integrals in (0) were evaluated 
taking s = 0. In the right panel of Fig.l we can-also see the behavior of the measured temperature that shows some 
qualitative agreement with experimental resultscl. 

Now, as we know the asymptotic behavior of the correlation and response, we can compute explicitly the large time 
evolution of the difference between the effective temperature and the measured one. A saddle-point calculation of (Q) 
shows that for t — > oo: 



f3 Y (t) 2E Y (t)f Y (t) 2 log(^L) \ 2 log(^ 



(15) 

In this calculation we made the assumption that due to the smallness of the relaxation time of the thermometer, 
the main contribution in the integral term of (0) is just when u = t. An important remark is that in this limit the 
difference does not depend on r provided that it is small enough (r/t << 1). This is clearly-seen in the inset of the 
left panel in Fig.l. If we consider only the second term in the r.h.s. of (||) (as proposed intl) we also find that the 
measured temperature converges to the effective one but the logarithmic asymptotic behavior in (|l5| ) is different. In 

g _« ( t) agy(t) 

fact the result for this case is " Z-} = F tt t- So, in order to obtain the correct logarithmic corrections to the 
asymptotic long-time behavior ( — > 0) we must consider the two terms in the expression for the current (|^). 
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FIG. 1. On the left there is a plot of the difference between the measured temperature and the effective one in function of 
the rescaled time (in Monte Carlo step unities) of measure. The thermometers have characteristic time-scales in the range: 
10 4 - 10 11 . The lowest curve corresponds to the shorter time-scale and the highest to the largest one. In the inset we have 
shown the asymptotic behavior for thermometers in the range: r = 10 3 - 10 4 being the upper line the analytical result (|l5|). 
On the right there is a plot of the measured temperature as function of the rescaled time for the same thermometers. The 
highest curve corresponds to the quickest thermometer. 



V. FOURIER'S LAW 



We have seen that the off-equilibrium current we have defined is a quantity which can be used to measure the 
effective temperature for systems which are off-equilibrium. Now we can think about the heat transfer itself. In 
standard non-equilibrium thermodynamics there is a linear relationship between the temperature gradient and the 
subsequent current establishecL-3. This is the well-known Fourier's law: 

J = Lqq\7(^) = -k\7T , (16) 

where Lqq is the Onsager coefficient and k is the thermal conductivity. For the equilibrium case we can compute the 
heat transfer (||) using (|l2|): 

limJ(M^4 /^y £y (fe-fr) , (17) 
t'—>t i\y 

Using the asymptotic expansions for /y(t) and Ey{t) (for a recent review seeEl) one can see that in the long-time 
limit, the non-equilibrium case satisfies Fourier's law: 

limJ(M')«e 2 - 1 l^ { p x .p Y{t)) . (is) 

(log(^)+5Mlog(^))) 2ilo ^(^) 

So the Onsager coefficient decays as Lqq ~ l/(f log(f)). As in the long-time limit the effective temperature (/3y(i)) 
decays like the energy, the thermal conductivity decays as n ~ log(t)/t. This simple example extends the Fourier's law 
to aging systems. So this is an extensiop-.of the linear relations between fluxes and forces obtained in the framework 
of the non-equilibrium thermodynamics!!-!], but with time-dependent Onsager coefficients. .— ■ 

Now, we can consider that system X is also a system of Monte Carlo harmonic oscillatoral3. In this case the 
Fourier's law reads (in the framework of LRT): 
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li m J(M0 = 8e ^^ /^ (fa-fr) . (19) 
t'-»t K X K Y Jx + JY 

In a recent papeiO we have shown that the effective temperatures and energies of two oscillator systems coupled 
like in (|J) and evolving under a sequential dynamics never equalize. So the heat transfer is extremely small and 
clearly insufficient te-thermalize degrees of freedom at different effective temperatures. The expression (|l9]) justifies 
the results found ir£3. Using the asymptotic expansions for fy{t) and Eylt) we can see that the Onsager coefficient 
decays as Lqq ~ l/(t log 2 (£)), and the conductivity as k ~ l/t. Following^ suppose that we couple the two oscillator 
systems at time t. If they are in the aging state evolving at zero temperature, then they have different effective 
temperatures given by A/3 = fix ~ Py ~ a l°g(i) with a a constant. The heat Q transferred between t and t + T is 
given by, 

V } J t x\og{x) SV log(i) > K ' 

For the two systems to equalize their temperatures the net heat transfer between t and t + T must be of the 
order Q ~ A/log(£) with A a finite constant. It can be shown that there are two regimes for Q depending whether 
T/t <C 1 or T/t ^> 1. If T/t <C 1 then Q ~ aT/(flog(t)) and the total heat transferred cannot equalize the effective 
temperatures. In the other regime the heat transfer can equalize the two temperatures only if T oc t. This implies 
that the time needed (in the aging regime) for compensate the difference between the two effective temperatures is 
of the order of the age of the system, then the heat transfer is unable to equalize the energies of the two systems in a 
finite time-scale so the zero-th law is not valid in the aging regime. 



VI. CONCLUSIONS 



We can now summarize our main results. We have established the proper way to measure the effective temperature 
of a glass. It is well known that the effective temperature based on the violation of FDT is a good quantity in 
the thermodynamic sense, so in the aging regime we can consider that the glass is in quasi-equilibrium at this 
effective temperature--. We have found that the correct procedure is to couple the thermometer weakly and measure 
the temperature of the thermometer which makes the heat transfer to vanish at a certain time t . In general, 
the measurement must be done with thermometers which have small characteristic time-scales compared to the 
decorrelation time of the aging system (in most cases this corresponds to the waiting time). 

To illustrate this, we have used a simple solvable model for a glass and we have computed the heat transfer established 
when we couple it with a generic thermometer. In this case the equations (||) and jlo|) were confirmed and also we 
have computed how the measured temperature converges to the effective one. Note that the model we have used is 
purely relaxational and only implies one large time-scale. Now, we know that for measuring the effective temperature 
we need a thermometer with a very short time-scale in comparison with the relaxational time-scale of the system. 
For the case of structural glasses is more complicated because in these systems there are two time-scales involved. 
The autocorrelation function reveals a secondary relaxation (/3-relaxation) and a primary one (a-relaxation) . For the 
secondary relaxation the system behaves as a system which is in equilibrium and FDT holds, whereas for the primary 
one, the relaxation of the system is extremely slow and FDT does not hold. Therefore, the effective temperature we 
have defined describes the violation of FDT in the long-time sector. If we are interested in the measurement of this 
temperature, our results suggest that we have to choose a thermometer with a time-scale shorter than the a-timescale 
but larger than the /3-timescale in order to avoid the effects of the initial /3-relaxation which is associated with a 
different temperature from the effective one. Otherwise, we would have a measurement where the two temperatures 
involved in the relaxation could mix. Along these lines it could_.be very instructive to investigate other models with 
two different time-scales (for instance the Backgammon modell__j). 

Employing our knowledge of the model of harmonic oscillators we have been able to extend Fourier's law to this 
glassy model by considering a time-dependent Onsager's coefficient; and we have seen how the thermal conductivity 
decays. It is important to note that the decay of the conductivity (in the case of two glasses coupled) is faster than 
the decay of the energy (i.e the decay of the effective temperature). Physically this means that, as time goes on, 
although the effective temperatures of the systems are different the conductivity is not high enough and they cannot 
thermalize. This seems to be the reason for the existence of different temperatures in real glasses, in which there are 
fast degrees of frjeedom and slow ones at different effective temperatures which never equalize due to the smallness of 
the conductivity!--.. 

It would be very interesting to address some otter non-equilibrium features such as the long-range spatial correlation 
functions for stationary non-equilibrium liquidsE-a for the model we have considered here and glassy systems in general. 
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